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Abstract 

We show that Zeeman excitations of ultracold Dy atoms trapped in an optical lattice can be used 
to engineer extended Hubbard models with tunable inter-site and particle number-non-conserving 
interactions. We show that the ratio of the hopping amplitude and inter-site interactions in these 
lattice models can be tuned in a wide range by transferring the atoms to different Zeeman states. 
We propose to use the resulting controllable models for the study of the effects of direct particle 
interactions and particle number-non-conserving terms on Anderson localization. 
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I. INTRODUCTION 


There is currently growing interest in engineering lattice Hamiltonians with ultracolcl 
atoms and molecules [T]. Of particular interest are extended Hubbard models, which include 
interactions between particles in different lattice sites. Such models exhibit rich physics and 
have been used to explain the role of long-range interactions in the context of superfluid 
- Mott insulator transitions [2], antiferromagnetism [3 0], high-Tc superconductivity [5], 
twisted superfluidity |6j, supersolids [7], self-trapping of bipolarons [8j. Extended Hubbard 
models are very difficult to solve numerically, especially for two- and three-dimensional 
lattices. Hence, the need to build experiments, where a many-body quantum system is 
described by an extended Hubbard model, whose parameters (in particular, the ratio of the 
hopping amplitude and the inter-site interaction energy) can be tuned by varying external 
fields, and where the particle densities can be imaged preferably with single site resolution. 
Tuning the parameters of the model, one could use such experiments to map out the phase 
diagrams. 

There are many proposals for realizing lattice models, including extended Hubbard mod¬ 
els [SHU], with ultracold atoms or molecules trapped in optical lattices. However, if ultracold 
atoms or molecules are used as probe particles of such models, the inter-site interactions are 
usually very weak. Therefore, the measurements of the phase diagrams require extremely 
low temperatures and extremely long coherence times, which are often difficult to achieve in 
current experiments. A more promising approach is to trap ultracold molecules in an optical 
lattice in a Mott insulator phase (with one molecule per site) and use rotational excitations 
of trapped molecules as probe particles of lattice models PEH23]. Such excitations can be 
transferred between molecules in different sites due to dipole - dipole interactions. The dy¬ 
namics of the excitations as well as their interactions can be controlled by external dc electric 
and/or microwave fields, leading to lattice models with tunable parameters. Experiments 
using excitations as probe particles of lattice models can tolerate much higher temperatures 
of atomic or molecular motion. However, it is currently not possible to create an optical 
lattice filled uniformly with molecules. On the other hand, ultracold atoms can be trapped 
in optical lattices with nearly uniform filling [21 2Tj. Thus, it would be desirable to engineer 
extended Hubbard models with internal excitations of atoms (instead of molecules) trapped 
in a Mott insulator phase. 


2 


A series of experiments have recently demonstrated the cooling of highly magnetic Cr 
[25] . Dy [25, [27], and Er [25J [2HJ atoms to quantum degeneracy. Such atoms interact via 
long-range magnetic dipole interactions and one can envision engineering the same lattice 
models with magnetic atoms as with ultracold molecules. However, the internal level struc¬ 
ture of magnetic atoms is more complex than the rotational structure of molecules and the 
nature of magnetic dipole interactions is different from that of electric dipole interactions. 
Motivated by the experiments on magnetic atoms and the work with ultracold molecules, we 
explore here the possibility of engineering extended Hubbard models with internal Zeeman 
excitations of ultracold magnetic atoms, such as Dy, trapped in a Mott insulator phase. 
Exploiting the unique nature of magnetic dipole interactions, we show that, for Zeeman 
excitations, the ratio of the hopping amplitude and inter-site interaction energy in the re¬ 
sulting lattice models can be tuned in a wide range by transferring the atoms to different 
Zeeman states. We discuss the advantages of using Zeeman excitations of magnetic atoms 
over rotational excitations of ultracold molecules. In particular, we show that the hopping 
of the Zeeman excitations in the lattice is insensitive to the magnitude of the magnetic field, 
which makes the coherent dynamics of excitations robust to field fluctuations. We show that 
Zeeman excitations in a diluted lattice of Dy atoms undergo Anderson localization over time 
scales less than one second and propose the models derived here for the study of the role of 
interactions and particle number fluctuations on Anderson localization. 
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II. LATTICE HAMILTONIAN WITH ZEEMAN EXCITATIONS 


We consider an ensemble of open-shell atoms with non-zero electron spin ( S ) and orbital 
angular momentum ( L) trapped in an optical lattice in the presence of an external DC 
magnetic held. We assume that the atoms fill the lattice uniformly with one atom per 
lattice site and that the atoms are not allowed to tunnel between different lattice sites. 
Thus, the atoms are separated by a large distance (> 260 nm) equal to half the wavelength 
of the trapping held. At such separations, the dominant interaction between the atoms 
in sites % and j is the magnetic dipole - dipole interaction V iy For simplicity, we assume 
that the atoms are arranged in a one-dimensional array along the z-axis of the space-hxed 
coordinate frame. In this case, 


a 


V* = 


2 L 


Jf | r/j T dy ./j_ . 


2 Jl. 


Jj,z 


( 1 ) 


where J z and J± are the ^-component and the raising/lowering operators of the total angular 
momentum J — L + S, acting on the space of the eigenstates \JM) of J 2 and J z , and a is 
the fine structure constant. The full Hamiltonian of the many-atom system is 


n = { AL * • + M-fci + 2 ^) ■ B } + (2) 

i i j^i 

where A is the constant of the spin-orbit interaction, //b is the Bohr magneton and B is the 
vector of an external magnetic held. 

We assume that all atoms are initially prepared in the Zeeman state | g) and a small 
number of atoms is then transferred to another Zeeman state |e). Note that the state \e) 
can be lower or higher in energy than the state | g). Following the approach described in 
Refs, plij (see also m), we derive the second-quantized Hamiltonian describing the Zeeman 
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transitions in this system: 


H ex = v g 


XX \ £e ' ~ £ 9 + X [( e il(^l^l e i)l^) - (9i\{gj\Vij\9i)\9j) > c f i,e'ke' (3) 
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7 l,e ,7 *.e" 7 i/' 


(9) 


where 




( 10 ) 


i j^i 


N is the number of atoms, e g and £ e / are the energies of the atomic states | g) and |e'), and 


Vi 9 / = {9i\{9j\Vij\9i)\9j)- 


( 11 ) 


Here, we assume that the Zeeman states ^ g and use the operators c\ e , and 

Cj, e ' defined by c\ e ,\gj) = %|e') and Cj >e '|e'-) = Sij\gj). For the purposes of this work, it is 
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convenient to rewrite this complex Hamiltonian as 


H ex — V g + y~~^(A£ eg + di)c\&i + ^ ' tijC^Ci + (12) 

* i 

^ X! U ij4 C i C j C J + ( 13 ) 

* 

\J2Y1 ( £ i a i + ^j) + ZJ s v (4 + c ') + S (4 + Ci) c )cj (14) 

* i j^i i i/* 

+U(e'^e,g) (15) 


where the operators c] and c t are defined by c]\gj) = Sij\ej) and c;|e g ) = 5ij\gj), Ae eg is the 
energy separation between the states \e) and | g), and the parameters of the Hamiltonian are 


dj - YAj, (16) 

d» = {Vff - V)f } , (17) 

% = Vff + V)f - 2Vff , (18) 

V‘f = V)f = (9 i |felhlft>K) (19) 

V“ = (e i |fe|V) j |e i )|e i ) (20) 

tij (fl l j| | |6j) |^j) (21) 

S H = ( e i\(9j\Vij\9i)\9j) (22) 


and 


P*i = ( e il(&l^l e i}l e l> - ( e i\(9j\Vij\9i)\9j)- 


(23) 


The terms (12), (13) and (14) are a part of the full Hamiltonian that describes the 
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Zeeman transitions only within the four-state subspace \a)\b) with both |a) and | 6 ) being 
either | g) or \e). If the energy gap for the | g) —y \e) transition were far detuned from all 
other energy gaps in the Zeeman level spectrum, it would be sufficient to consider the part 


of the Hamiltonian given by Eqs. (12), (13) and (14). It is important to note that for highly 
magnetic atoms it may be necessary to consider Zeeman states outside of this subspace. 
Figure la shows that the Zeeman states of a Dy atom in the ground electronic state form 
a ladder of nearly equidistant levels at weak magnetic fields. This pattern of energy levels 
is characteristic of highly magnetic atoms with zero or negligible hyperhne structure. This 
pattern of energy levels allows for transitions to states outside of the subspace spanned by 
\g) or |e). For example, two atoms in the | g) state may interact to produce two Zeeman 
states with energies just above and just below that of | g). Such interactions are induced 
by the matrix elements in Eq. ([7]). The full Hamiltonian must also include the terms that 
describe the interactions of two atoms in states e', e" ^ g to produce atoms in other states 
/', f" 7 ^ g, e. Since the majority of atoms are in a particular state | g), we assume that such 
interactions are unlikely and neglect them. 

Various lattice models can be engineered by controlling the magnitude of the different 
matrix elements of the magnetic dipole interaction entering Eqs. 


III. ENGINEERING LATTICE MODELS 

In this section we show (i) how to simplify the lattice Hamiltonian presented in Section II 
by applying magnetic fields; and (ii) how to tune the relative magnitudes of the parameters 
of the resulting lattice models by transferring atoms into different states. We illustrate the 
tunable range of the parameters by calculating the model parameters for the specific example 
of Dy atoms in an optical lattice. 


t — V model 


Eqs. (12) and (13) represents a t — V model (32], an extended single band, Hubbard model 


for hard-core bosons da lani eh. This model can be studied with the Zeeman excitations if 


the effect of the terms (14) and (15) are suppressed. As we show below, this can be achieved 


by applying a finite magnetic held and introducing a small admixture of different M -states 
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into the eigenstates | JM). 


Eqs. ([3]) - ([9]) and (12) ~ (14) can be separated into terms that conserve the number of 
excitations (Eqs. [3] —15, 12 and 13) as well as particle number-non-conserving terms (Eqs. [b] 
— [ 9 ] and 14). If the Zeeman states form a ladder of equidistant states, the particle number- 
non-conserving terms can be further separated into energy-conserving (some terms in Eq. [7]) 
and energy-non-conserving terms (Eqs. 3 BE The effect of the energy-non-conserving 
terms can be eliminated by applying a finite magnetic field such that the energy difference 
between the Zeeman levels is significantly larger than the magnitude of the matrix elements 
appearing in Eqs. (g-g) and ( 0 . 

In order to eliminate the effect of all terms in Eq. (15), it is necessary to make the energy 
gap for the | g) —» |e) transition unique, i.e. different from the energy gaps in the Zeeman 
spectrum just below and just above the states | g) and |e). This can be achieved by applying 
a magnetic field strong enough to shift the Zeeman levels due to couplings between different 
total angular momentum states. As illustrated in the lower panel of Figure 1, these couplings 
introduce a differential in the energy gaps between different Zeeman states. To illustrate 
this, we plot in Figure lb the of the energy gaps between the states correlating with the 
states | J — 8, M — —1) and | J = 8 ,M — 0); states \J — 8, M — 0) and | J — 8, M — +1) 
and states | J = 8, M = +1) and \J — 8, M — +2), as functions of B 0 . As Figure lb shows, 
the magnetic field with B 0 & 200 — 300 G produces the differential of the energy gaps equal 
to the matrix elements t^i + 1 for Dy atoms on an optical lattice with a = 266 nm. At fields 
with B 0 > 300 G, the difference in the energy gaps becomes much larger than any of the 


matrix elements in Eq. (15) so the Hamiltonian (12) - (15) reduces to the t — V model. 

The parameters of the t — V model can be tuned by transferring atoms into different 
Zeeman states. If the | g) and |e) states are the Zeeman states | g) = | JM) and |e) = \JM'), 
the matrix elements (11) and ([21]) of the operator ([l]) can be written as follows: 


and 


= y ge - V" 


— ( M 2 - M'M) 
V ij 


(24) 


Q> . ... 

tij = bTT M’,M-1 + a -^+^M',M-lKl',M+l] 


( 25 ) 


















with 


<4 = W + 1) - M(M ± 1)] 1/2 
V ± = [J(J + 1) - M'(M' ± 1)] 1/2 


(26) 

(27) 


The interaction between the Zeeman excitations (18) can be written as 


V» = - [«? - V)f ) + (vg» - V*)} =~(M~ M'f 


(28) 


v 


These equations show that the diagonal matrix elements V 99 and and hence d l3 and 
Vij are non-zero, provided both M/0 and M' ^ 0. This is different from the case of the 
electric dipole - dipole interaction between molecules [35] . The electric dipole interaction 
must couple states of the opposite parity. Therefore, if | g) and |e) are the eigenstates of 
a molecular Hamiltonian in the absence of electric fields, the matrix elements dij and v t] 
of the electric dipole - dipole interaction vanish. These interactions can be induced in an 
ensemble of polar molecules by applying an external electric field that mixes the rotational 
states with different parity 


In contrast, the matrix elements of the magnetic dipole - dipole interaction (24) and (25) 


should not be expected to vary significantly with an external magnetic field. This will be 
illustrated and discussed in the following section, using the example of Dy atoms on an 


optical lattice. As follows from Eqs. (24) and (25), the relative weights of the two couplings 


can be tuned by choosing different Zeeman states | JM) as the | g) and \e) states. Notice, 
for example, that for the particular case of \g) being the state | J,M = 0), the magnitudes 
of , and consequently di, vanish. 


B. t — V model with Dy atoms 

We illustrate the range of controllability of the parameters of the t — V models using 
an example of Dy atoms in an optical lattice. The absolute magnitudes of d i3 , tij and Vi 3 
increase with J as the square of the magnetic moment. The ground electronic state of Dy 
is characterized by the total angular momentum J = 8 so Dy atoms have a large magnetic 
moment (10 Bohr magnetons) and a manifold of Zeeman states displayed in Figure la. The 
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Zeeman structure of Dy allows for the possibility of using the state |M = 0) as the | g) state, 
leading to the value dij = 0. 

If the states for the Zeeman excitations in an ensemble of Dy atoms are chosen to be 


well-defined angular momentum states \g ) = | JM) and \e) = \JM'), Eq. (25) shows that 
tij = 0 unless |M — M'\ = 1. Eq. (28) shows that the interaction Vij is oc (M — M') 2 so it 
is independent of M and M'. if \M — M'\ = 1. However, the parameter t l3 is sensitive to 
the magnitudes of M and M'. This is illustrated in the upper panel of Figure 2. The ratio 
t ^/can thus be tuned by transferring atoms into the Zeeman states with different M, as 
illustrated in the lower panel of Figure 2. Notice that the ratio Uj/vij is always negative, 
which means that the interactions between the excitations are always effectively attractive. 


The largest magnitude of the ratio tij/vij ~ —18 can be achieved when the atoms are 
prepared in the Zeeman state with M — 0 and excited to the Zeeman state with M = +1, 
while the smallest magnitude of the ratio t l3 /v i3 & —4 can be achieved by preparing the 
atoms in the maximally stretched state \J — 8, M = —8) or \ J — 8, M — +8). 

As illustrated in Figure 4 the absolute magnitude of v l3 can be tuned if the atoms are 
prepared in coherent superpositions of states with different M. Consider for example the 
superpositions | g) = a\JM) + f3\J,M + 5) and \e) = a!\JM') + /3'\J,M' + S'). For the 
parameter t l3 to be non-zero, either | M — M'\ or | M — M' + 5 — T| must be 1. However, 
there is no such restriction on the matrix elements determining the magnitude of ty,-. As 


follows from Eq. (28), the magnitude of is expected to increase with increasing the 
difference between the angular momentum projections of the states participating in the 
excitation. This is graphically illustrated in Figure 3, showing that the magnitude of v t j can 
reach 600 Hz, if M — M' = 16. This suggests that the ratio jv l3 can be tuned by preparing 
the atoms in the coherent superpositions of the following kind: | g) = a\JM) + /3| J, M + S) 
and \e) = a'\JM + 1 ) + /3'\J, M + 5'). The parameters t l3 and v l3 for these states are both 
non-zero and the magnitude of j can be modified by varying the value of |<5 — 5'|. 

The interaction of atoms with a magnetic field couples states with different total angular 
momenta J, which may - in principle - modify the atomic states | g) and |e), and, conse¬ 
quently, the lattice model parameters. It is important to examine the effect of an external 
magnetic field on the lattice model parameters. To do this we diagonalized the full Hamil¬ 
tonian of the Dy atom in a magnetic field B = B 0 (O.lx + z) and used the eigenstates to 


evaluate the model parameters in Eqs. (12) - (13). Since the states of different J in the Dy 
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atom are separated by large energy gaps (> 1000 cm -1 ) due to the spin-orbit interaction, 
the eigenstates of Dy in a magnetic held are nearly identical to the angular momentum 
states | JM). Figure 5 shows the nearest-neighbour coupling parameters and v iji+ i for 
a one-dimensional array of Dy atoms on an optical lattice with the lattice site separation 
a = 266 nm computed for two pairs of Zeeman states at different magnetic fields. The 
results shown in Figure 4 illustrate that the Hamiltonian parameters do not change with the 
magnetic held in the interval of held strengths between zero and 5000 G. This is important 
because it shows that the magnetic held can be used to separate the Zeeman states in order 
to create isolated two-level systems or tuned to the limit of vanishing held where the terms 


in Eq. (14) become important, without affecting the parameters of excitation interactions. 


C. Particle number-non-conserving interactions 


In the limit of weak magnetic helds, as Ae eg —> 0, the energy separation between different 


particle number states of the model (12) decreases to the minimum of di. As follows from 


Eq. (24), this parameter can be eliminated if the ground state | g) is chosen to be \J,M = 0). 


At weak magnetic helds, the particle number-non-conserving terms (14) as well as the terms 


in Eq. (15) must be included in the Hamiltonian. Number non-conserving interactions may 


mediate effective long-range hopping (for example, a particle can move in a lattice by virtual 
transitions to the three-particle subspace and back). As such, these interactions may have 
non-trivial effects on the dynamics of quantum walks and localization in disordered lattices. 
Such interactions arise in the context of excitons in molecular crystals [37]. However, they are 
usually negligibly small and difficult to investigate. As shown below, number non-conserving 
interactions can be made significant in the system considered here. 

We hrst note that if the array of atoms is arranged along the magnetic held direction, 
the matrix elements of the operator (JTj) that determine the parameters s t j and pij in Eq. 


(14) vanish. This simplihes the resulting lattice models to include only the hrst of the 


particle number-non-conserving terms in Eq. (14). If desired, the terms s t j and pij can be 
tuned to Unite values if the magnetic held direction is changed or the atoms are prepared 
in coherent superpositions of different M-states. For example, if | g) = \J,M = 0) and 
\e) = a\J = 8, M = 0) + (3\J — 8, M = 1), all of tij, Sij and p t j become non-zero. Here, 
we assume that the magnetic held is directed along the atomic array and that s tJ = 0 and 
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Vij = 0 . 

Care must be taken when considering the limit A e eg —y 0. In this limit, multiple Zeeman 
states become degenerate and it may be necessary to consider interband couplings deter¬ 


mined by Eq. (15). This may be useful if complicated models, including multiple excitations 
of different kind, are desired. Note, however, that if | g) and |e) are states with well-defined 
M and M' , a two-atom state \M)\M') can only be coupled to the same state, the state 
\M')\M) or a state | M ± 1)| M' =|= 1). The matrix elements of the dipole - dipole interaction 
(M, M' | V t j | M ± 1, M' =p 1) change the number and type of excitations in the atomic ensem¬ 
ble. These processes can be eliminated if the state | g) is chosen to be | J,M = ±J). In this 
case, the effective lattice model describing the dynamics of | g) —> \e) excitations is 


H ex Vg T y ^ (A £ e g T (l ;) c,j Cj T ^ 'y ^ UjCjCi T 


* 


v v c h c h + + d A 


* 3 # 


i j^i 


(29) 


It is important to note that this model is valid as long as Ae eg (which is determined by the 
magnitude of the magnetic held) is significantly larger than tij. In this limit, the effect of the 
number-non-conserving terms is perturbative, i.e. a single excitation remains predominantly 
in the single-particle subspace, undergoing virtual transitions to the three-particle subspace. 


If the energy gap A e eg is so small that the interactions (29) lead to the creation of multiple 


particles, other terms in Eq. ([?]) must be included, making the Hamiltonian more complex. 

If the effects of the interactions are to be removed, one can choose the states | g) = 

| J,M = 0) and \e) = |J, M — 1). In this case, \tij\ \vij\ (see Figure 4). However, the 
lattice model for these excitations is also affected by terms in Eq. (J7|), which lead to leaking 
of the |e)-state populations to other Zeeman states of higher energy. These terms lead to 
the spontaneous creation of atoms in Zeeman states above and below the energy of the state 
with M — 0, as well as the inverse process. The Zeeman state populations must eventually 
return to states e and g, as the total number of the Zeeman states is finite and small. These 
terms thus serve as an additional source of particle number-non-conserving interactions that 
generate atoms in state e. 
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D. Anderson localization of Zeeman excitations 


Until now, we assumed that the atoms populate the optical lattice uniformly. If the 
lattice is populated partially (which is more often the case in experiments), the empty lattice 
sites serve as impurities that can scatter the Zeeman excitations. Since the distribution of 
empty sites is random, the Zeeman excitations thus propagate in a randomly diluted lattice. 
Tuning the models as described above suggests an interesting opportunity to explore the role 
of direct particle interactions and number non-conserving forces on Anderson localization 
in disordered lattices [311,1321 • I n addition, the ability to design optical lattices with various 
dimensionalities and geometries can be used to verify the scaling hypothesis of Anderson 
localization ra as well as Anderson localization of particles with long-range hopping in 
various geometries m Here, we explore if the parameters of the models based on Zeeman 
excitations of Dy are significant enough to allow Anderson localization over experimentally 
feasible time- and length-scales. 

We consider an isolated Zeeman excitation in a one-dimensional array of 1000 Dy atoms 
trapped in an optical lattice with a = 266 nm containing 20 % of empty lattice sites. We 
use the parameters correpsonding to the | J = 8,M = 0) — y \J = 8 ,M = +1) excitation 
and compute the dynamics of quantum walk for the Zeeman excitation placed at t — 0 on a 
single atom in the middle of the lattice. The wave packet of the excitations is propagated by 
computing the time-evolution operator, as described in detail in Ref. [32] • The results of each 
dynamical propagation are averaged over 100 disorder realizations (random distributions of 
empty lattice sites). 

The results shown in Figure 6 illustrate that the Zeeman excitation forms an exponentially 
localized spatial distribution within one second. The width of the distribution characterized 
as the length L containing 90 % of the excitation probability exhibits a short-time oscillation 
which is likely an effect of coherent back scattering and approaches the value of ~ 20 lattice 
sites in the limit of long time. These results can be directly mapped onto the results 
describing Anderson localization for rotational excitations in an ensemble of polar molecules 
ra and the electronic excitations in an ensemble of Rydberg atoms pf3 j. 
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IV. CONCLUSION 


In this work, we consider Zeeman excitations in an ensemble of highly magnetic atoms 
(such as Dy) trapped in an optical lattice, with one atom per lattice site. The Zeeman 
excitations can travel in the lattice due to energy transfer between the atoms. The most 
important results of this work can be summarized as follows: 

• We show that superpositions of the Zeeman excitations can be used to simulate the 
t — V model (the single-band, extended Bose-Hubbard model for hard-core bosons). 
The parameters of the model (most importantly, the ratio of the hopping amplitude 
and the inter-site interaction energy) can be tuned by preparing the atoms in different 
Zeeman states. For an ensemble of Dy atoms on an optical lattice with a = 266 nm, 
we show that the inter-site interaction can be engineered to be as large as 600 Hz. 

• We illustrate that the parameters of the model (hopping amplitudes and inter-site in¬ 
teractions) are insensitive to the magnetic held. This has two significant consequences. 
First, an external magnetic held can be used to uncouple the electron degrees of free¬ 
dom from nuclear spins, thereby removing complications associated with the hyperhne 
structure of atoms and the degeneracies of the Zeeman states. Second, an external 
magnetic held can be used to separate the Zeeman states, leading to suppression of 
energy- and particle number-non-conserving terms. 

• We show that the same Hamiltonian can be used to simulate a lattice model with 
significant cjcj terms, leading to particle number interactions. These interactions 
mediate effective interactions modifying the hopping of particles and can be used to 
produce entangled pairs 031. 

• Since the lattice with randomly distributed empty sites leads to a quantum percolation 
model for the Zeeman excitations, we propose to apply the models derived here for 
the study of Anderson localization induced by off-diagonal disorder. In particular, 
our results suggest the possibility of studying the role of inter-site interactions and 
particle number huctuations on quantum localization in diluted lattices. We show that 
for an optical lattice with a = 266 nm partially populated with Dy atoms, Anderson 
localization of excitations placed on individual atoms occurs over timescales less than 
a second. 
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FIG. 1: Upper panel: Zeeman levels of a Dy( 5 /) atom in the lowest-energy spin-orbit state char¬ 
acterized by J = 8 in a magnetic field B = Bqz. Lower panel: the solid curve - difference of the 
energy gaps (em =2 — £m= 1 ) — (em =1 — £m= 0); the dot-dashed curve - difference of the energy gaps 
(em= 2 — £m= 1 ) — (em =0 — £m=-i)- The horizontal dashed line shows the magnitude of the matrix 
element in Eq. (21) for Dy atoms with Ig) = \ J = 8, M = 0) and |e) = \ J = 8, M = 1) in an 
optical lattice with a = 266 nm. 
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FIG. 2: The magnitudes of the coupling constants fjj (upper panel) and the ratio tij/vij (lower 
panel) with j = i ± 1 for the Zeeman states of Dy corresponding to | g) | JM) and |e) | JM'). 

The calculations are for the magnetic field B = Bq (O.lx + z) with Bq = 100 G. The Zeeman states 
in this magnetic field retain 96% of the eigenstates of J 2 and J z . 
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FIG. 3: The magnitude of the coupling constant Vij with j = i ± 1 for the Zeeman states of 
Dy corresponding to | g) =f- \JM) and \e) => The calculations are for the magnetic field 

B = Bq (O.lx + z) with Bo = 100 G. The Zeeman states in this magnetic field retain 96% of the 
eigenstates of J 2 and J z . 
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FIG. 4: The magnetic field dependence of the quantities (squares) and i^j+i (circles) defined 
in Eqs. (24) and (25) for two different pairs of the Zeeman state of Dy(J = 8) atoms: the full 
symbols - the results for | g) = |J = 8, M = —8) and |e) = \J = 8, M = —7); the open symbols - 
the results for \g) = \J = 8,M = 0) and \e) = \J = 8, M = +1). The magnetic field is given by 
B = Bq (O.lx + z). The Zeeman states in such a magnetic field retain 96% of the eigenstates of 
J 2 and J z . 
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FIG. 5: The ratio for the Zeeman states of Dy corresponding to \g) = \ J = 8, M = —7) 

and |e) = a\J = 8, M = —8) + b\ J = 8, M'): full circles - a = 1/4; open circles - 

a = \/3/5, b = \/2/h. The calculations are for the magnetic field B = Bq (O.lii + z) with Bo = 100 
G. The Zeeman states in this magnetic field retain 96% of the eigenstates of J 2 and J z . 
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FIG. 6: Anderson localization of the | J = 8 ,M = 0) —> \ J = 8 ,M = +1) excitation in a one¬ 
dimensional array of Dy atoms on an optical lattice with a = 266 nm and 20 % of the lattice sites 
empty. The upper panel shows the probability distribution for the atoms in the corresponding site 
to be in the excited state at t = 2 seconds formed by a single excitation placed at t = 0 in the 
middle of a lattice with 1000 sites. The lower panel shows the width of the excitation probability 
distirbution as a function of time. 
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